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Abstract. We study a general stochastic rumour model in which an ignorant individual 
has a certain probability of becoming a stifler immediately upon hearing the rumour. We 
refer to this special kind of stifler as an uninterested individual. Our model also includes 
distinct rates for meetings between two spreaders in which both become stiflers or only one 
does, so that particular cases are the classical Daley-Kendall and Maki-Thompson models. 
We prove a Law of Large Numbers and a Central Limit Theorem for the proportions of 
those who ultimately remain ignorant and those who have heard the rumour but become 
uninterested in it. 



1. Introduction 

Various mathematical models for the propagation of a rumour within a population have 
been developed in the past decades. Two classical models were introduced by Daley and 
Kendall ^ and Maki and Thompson [9|. In the first model (which we shorten to [DK] 
model), a closed homogeneously mixing population of iV + 1 individuals is considered. Peo- 
ple are subdivided into three classes: ignorants (those not aware of the rumour), spreaders 
(who are spreading it), and stiflers (who know the rumour but have ceased communicating 
it after meeting somebody who has already heard it). We adhere to the usual notation, 
denoting the number of ignorants, spreaders and stiflers at time t by X{t), Y{t) and Z[t), 
respectively. Initially, X(0) = N, F(0) = 1 and Z(0) = 0, and X{t) + Y{t) + Z{t) = N + 1 
for all t. The process {(X(t), Y{t))}t>o is a continuous-time Markov chain with transitions 
and corresponding rates given by 

transition rate 

(-1,1) xr. 



(0,-2) 



'Y' 

(0,-1) Y{N + 1- X -Y). 

People interact by pairwise contacts, and the three possible transitions correspond to 
spreader- ignorant, spreader-spreader and spreader-stifier interactions. In the first case. 
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the spreader tells the rumour to an ignorant, who becomes a spreader. The two other 
transitions represent the transformation of the spreader(s) involved in the meeting into 
stifler(s), that is, the loss of interest in propagating the rumour derived from learning that 
it is already known by the other individual in the meeting. 

In the model formulated by Maki and Thompson (of (which we shorten to [MT] model), 
differently from the [DK] model, the rumour is spread by directed contact of the spreaders 
with other individuals, i.e., the initiator is distinguished from the recipient. In addition, 
when a spreader contacts another spreader, only the initiating one becomes a stifler. Thus, 
the continuous-time Markov chain {{X{t),Y{t))}t>Q evolves according to the following 
table 

transition rate 
(-1,1) XF, 
(0,-1) F(iV-X). 

We refer to Daley and Gani [2, Chapter 5] for an excellent account on the subject of rumour 
models. 

We introduce a general stochastic rumour model with the following characteristics: 

(i) Distinct rates for meetings between two spreaders in which both become stiflers or 
only one does. 

(ii) A new class of individuals, which we call uninterested. In traditional models, once 
an ignorant is told the rumour, only the transformation into a spreader is allowed. 
In our model, this individual has the choice (with a certain probability) of becoming 
uninterested, that is, of stifling right after hearing the rumour. 

We underline that a stifier is an ex-spreader who has lost interest in propagating the 
rumour, whereas an uninterested is an ex-ignorant whose lack of interest arose immediately 
upon hearing the rumour from a spreader. 

To define the model, we consider the notation introduced so far and denote by U {t) the 
number of uninterested individuals at time t. Initially, X{0) = N, f/(0) = 0, Y{0) = 1 and 
Z(0) = 0, and X{t) + U{t) + Y{t) + Z{t) = iV + 1 for all t. Let V{t) = (X(t), U{t), Y{t)). 
We suppose that {V{t)}t>Q is a continuous-time Markov chain with initial state (X, 0, 1) 
and 

transition rate 
(-1,0,1) X6XY, 
(-1,1,0) A(1-5)XF, 

(0,0,-2) Ar'^' 
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(0,0,-1) X92Y{Y -1) + XjY(N + 1-X -Y). 

The first two cases correspond to a spreader telling the rumour to an ignorant, who decides 
to become a spreader or an uninterested (that is, an immediate stifier) with respective 
probabilities 6 and 1 — 6. The interaction between two spreaders splits into two cases, 
in which both become stiflers or only one does. These cases are represented by the third 
transition and by the first part of the rate of the fourth transition. Finally, the second part 
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of this rate corresponds to a meeting between a spreader and a stifler or an uninterested 
individual. 

We define 6 = 6i + 62 — ^ and assume throughout the paper that 

(1.1) A > 0, 7 > 0, ^1 > 0, ^2 > 0, < 5 < 1 and < ^ < 1. 

Notice that both [DK] and [MT] models can be obtained by suitably choosing these con- 
stants. The family of Markov chains just defined includes the classical and other ru- 
mour models, some of which are presented in ^ showing the convenience of the adopted 
parametrization. 

To state our results, it is convenient to write down explicitly the dependence of V on A^, 
so we use the notation V'^^'^ = {X^^\U^^\Y^^^). Observing that the process eventually 
terminates (when there are no spreaders in the population), we define 

r(^) = inf{t:y(^)(t)=0}. 

In ^ we present a Weak Law of Large Numbers and a Central Limit Theorem for 
N^^ (X*^^)(r'^^)), [/^^^(r'-^-')), that is, the ultimate fractions of the originally ignorant in- 
dividuals who remained ignorant and who became uninterested. Proofs are presented in 

m 

As far as we know, limit theorems were rigorously proved only for the basic [DK] and 



[MT] models (see Pittel [Jjj, Sudbury [l2j and Watson [13[). In addition, the main tools 
applied for stochastic rumour processes have been the analysis of the embedded Markov 
chain, martingale arguments, diffusion approximations, generating functions and the study 
of analogue deterministic versions. A good survey of these methods can be found in Daley 
and Gani [2|. Briefly, the technique we use consists in defining a coupled process with 
the same transitions as the original process until absorption. This construction is suitably 
defined, in such a way that the new process is a density dependent stochastic model for 
which the theory presented in Ethier and Kurtz [3] can be applied. This approach allows 
us to establish limit theorems for more general models. To the best of our knowledge, 
this idea is used for the first time in the context of stochastic rumour models here and in 
Lebensztayn et al. jl], in which we study a generalization of the Maki-Thompson model 
with random stifling and general initial configuration. 



2. Main results 
We start off with a few definitions. 

Definition 2.1. For < < 1, consider the function / : [0, 1] — ^ M given by 

_ (7 + S9)x' - (7 + 6)9x - 7(1 - 9) 
9(1-9) 

We define x^o = Xoo(^, 7,^) as the unique root of f{x) = in the interval (0, 1). 
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To justify the existence and uniqueness of the root, notice that /(O) < 0, /(I) = and 
that / is unimodal, with a global maximum at the point 

-' + *''^'^e(o,i). 



7 + 5 

For 9 = 1/2, an explicit formula for be obtained, namely, 

a;oo(5,7,l/2)= 

In the cases 6 = and 6 = 1, x^o is defined similarly. 
Definition 2.2. Consider the functions /o and fi defined on (0, 1] by 
fo{x) = lim f{x) = (7 + (5)(1 - x) + 7logx, 

= lim /(x) = —7(1 — x) — {-y + 6) X logx. 

For each 6 G {0, 1}, we denote by Xoo{S,'y,6) the unique root of fg{x) = in (0, 1). 

Remark 2.3. In the last be written in terms of the Lambert W function, which 

is the multivalued inverse of the function x ^ xe^. Let Wq and W^i be the principal and 
the lower real branches of the Lambert W function, respectively; see Corless et al. [ij for 
more details. Defining /i = 1 + 5/7, we have that 

x^{6, 7, 0) = -h-^ Wo{-h e~^) and 

x^S,^,!) = -[hW^^i-e-'/''/h)Y\ 

We prove that, whatever the value of 6, the following inequality holds: 

(2.1) Xoo(5,7,^) < 

7 + 

For this, it is enough to show that / and fg evaluated at the point 7/(7 + 5) are strictly 
greater than zero. For 6 G {0, 1}, formula (12. ip is therefore an immediate consequence of 
a standard logarithm inequality: 

u — 1 

(2.2) < logM < M - 1 for all < M < 1. 

u 

For 6 G (0, 1), we use (12. 2p to prove that the function 6 1— )• (7/(7 + 56'))^/^ is increasing, 
whence (12. ip follows. 

We are now ready to state our main results. 

Theorem 2.4. Assume (II. ip and let x^o be given by Definition \2.1\ or \2.2 according as 
6 G (0, 1) or not. Define u^o = (1 ~ 5)(1 ~ a^oo)- Then, 

XW(rW) , f/W(rW) 
hm = Xao and lim = tn probability. 

N^oo N N^oo N 
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Theorem 2.5. Assume (11 .ip and define 



K = 39i + 202 - 47, A 



B 



'jSUc 



C = (7 + 5)2 (45^2 _ ^ 266)) + K7(7 + 5){^ + 6(29 - 1)) - 4^72(1 



C{l-x^) 



27 



2(20- l)(7 + 5^)^ 



«:5(27 + 5) - 27(5 - /s:(7 + 5)) log 



(7 + 6y 



N 



f/W( 



r 



Mr 



V 



D = 

Then, 
(2.3) 

where — )■ denotes convergence in distribution, and N2{0,T,) is the bivariate normal distri- 
bution with mean zero and covariance matrix 

S21 S22 



/XW(r(^)) 



^ A/'2(0,S) as ^ cx), 



whose elements are given by 

Sll = Xc 



(2.4) 



S12 
S22 



^21 



-Xoo) + A^D, 
-(l-5)Sn + A5, 



(1 - 5)2Sn + (1 - - Xoo) - 2AB). 



Of course, when 5 = 1, we can omit the second component of the vector in the left- 
hand side of (12. 3p . and we denote by (= Sn) the variance of the asymptotic mean zero 
Gaussian distribution. 

Example 2.6. Let p G [0, 1] and consider our model with the choice A = 5 = 7 = l, 6'i=p 
and 6'2 = 1 — p, so = 0. Thus, the limiting proportion of ignorants and the variance of 
the asymptotic normal distribution in the CLT are given respectively by 



Xoo(l, 1,0) 



iyo( -Se- 



er 



a;oo(l - a;oo) (1 - 2 + 2 px; 



1 



2 3^00)^ 



0.203188, and 



0.272736 + 0.0379364 p. 



We obtain the [MT] or [DK] model according to whether p equals or 1, showing that our 
theorems generalize classical results presented by Sudbury IJ] and Watson 13|. 



Example 2.1 . In an interesting confessional essay, Hayes [5j describes a mistake committed 
when he attempted to simulate the basic [DK] model. The author actually simulated the 



6 



ELCIO LEBENSZTAYN ET AL. 



[MT] model, with the difference that, when two spreaders meet, both become stiffers. This 
model is obtained by choosing A = 5 = 7 = l, 0i = 2 and 62 = 0, in which case = 1, 

Xoo = a^oo(l, ^' ^) = ~ 2 ^(_Vv2/2) ^ 0-284668, and 

a' = ^00(1- Xoo) (l-3xoo + 3a;g,) ^ q_427204. 
(1-2X00)2 

This clarifies the numerical value of the proportion of ignorants remaining in the population 
that Hayes obtained in his simulations. 

Example 2.8. Let a,p,q & (0, 1]. We describe a new variant of the [DK] model, which we 
call the (q;,p, g)-[DK] model. Suppose that, independently for each pairwise meeting and 
each individual, 

(a) A spreader involved in a meeting decides to tell the rumour with probability p. 

(b) Once such a decision is made, any spreader in a meeting with somebody informed 
of the rumour has probability a of becoming a stifler. 

(c) Upon hearing the rumour, an ignorant becomes a spreader or an uninterested in- 
dividual with respective probabihties q and 1 — g. 

This model with q = 1 (no uninterested individuals) was introduced by Daley and 
Kendall [sj, who studied its deterministic analogue. This analysis is also presented in 
Daley and Gani [1, Section 5.2]. The basic [DK] model corresponds to a = p = q = 1. 

Observe that, for the (a,p, g)-[DK] model, the continuous-time Markov chain {V(t)}t>o 
evolves according to 

transition rate 

(-1,0,1) pqXY, 

(-1,1,0) pil-q)XY, 

(0,0,-2) «M2-P)'^ 



2 

(0, 0, -1) a{l - a)p{2 - p) Y{Y - 1) + apY{N + 1-X-Y). 



Therefore, Theorems 12.41 and 12.51 yield novel bivariate limit theorems for this model, which 
are obtained by making the following substitutions: 

\ = p, 5 = q, 61 = a^{2 — p), 62 = a(l — a){2 — p), 'j = a and 9 = a(l — p). 

Here are some important cases: 

(1) For the {a,p, g)-[DK] model with p < 1, the asymptotic proportion of ignorants Xoo 
is the unique root of 

^ (1 + qjl - p))x"(i-P) -(a + q){l-p)x-l + a{l - p) ^ 

^^""^ {l-p){l-a{l-p)) 

in the interval (0, 1). The limiting fraction of uninterested individuals is u^o = (1 — g)(l — 
3^00)- When 5=1, the value Xoo coincides with that obtained for the deterministic analogue 
of the model in Daley and Gani [2, Equation (5.2.8)]. 
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(2) The (l,l,g)-[DK] model is the basic [DK] model with the additional rule that an 
ignorant is allowed not to have interest in spreading the rumour. For this model, the 
limiting proportions of ignorant and uninterested individuals are expressed respectively as 

, ^. Xoo = a;oo(g, 1,0) = — ; with/i = l + g and 

(2.5) h 

Moo = (1 - g)(l - a;oo)- 
The entries of the covariance matrix S given in (12.40 simplify to 

a;oo(l - Xoo) (2 - (3 + g2)a;^ + (1 + g)2a;2^) 



S 



2(l-(l + g)xoo)2 
Xoo (-2(1 - g) + (1 - g)(3 + q)x^ - (1 + qfxi;) 



12 



2(1- (l + g)a; 



^2 



oo I 



y ^ Uoo (2g + 2(1 - 5g)xoo + (-3 + 9q + 3q^ - q^)xl + (1 - g)(l + qfxl) 

2(l-(l + g)xoo)2 

When q = 1, these formulae reduce to the well-known results for the basic [DK] model. 
(3) For the {a, 1, 1)-[DK] model, we have that 

Xoo = 3;oo(l, 0) = ; with h = 1 -\ — , and 

h a 



2 



a 



Xoo(l - Xoo) (2«^ + [2(1 - a) - a{l + a^] Xoo + a(l + a)^a;^) 

2(a - (1 + «)xoo)2 ' 



These formulae agree with those presented in Exercise 5.7 of Daley and Gani ^] (in which 
the reader is asked to obtain by making use of Kendall's Principle of Diffusion of 
Arbitrary Constants). 

Example 2.9. We define the {a,p,q) version of the [MT] model in a similar way, with the 
rules (a) and (c) given in Example 12.81 and 

(b') Once a spreader decides to tell the rumour in a directed contact with somebody 
already informed, only this spreader chooses with probability a to become a stifler 

instead of (b) in Example 12.81 holding independently for each directed contact between 
two individuals. The basic [MT] model has a = p = q = 1. 
Since the {a,p,q)-[MT] model is obtained by the choice 

\ = p, 6 = q, 9i = 0, 02 = 'J = a and 9 = 0, 

we conclude that the limiting proportion of ignorants is given by 

Xoo = Xoo{q, ct, 0) = — — ^ with h = 1 + —. 

h a 

Two particular cases are: 
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(1) For the (1, l,g)-[MT] model, the values of x^o and the same as those of the 

1^1, 1, g)-[DK] model, given by (12. 5p . However, 

Xoo{^ - Xoo) (1 - (1 + g^)xoo) 



^22 



(l-(l + g)xoo)^ 



2 



(l-(l + g)xoo)2' 
Moo + (1 - 5g)xoo + (-1 + 4g + 



2^^2 ^ 

oo/ 



(l-(l + g)xoo)2 

(2) The (a, 1, 1)-[MT] model corresponds to the basic [MT] model in which the number 
of contacts with an informed individual by each spreader waiting to become a stifler has 
geometric distribution with parameter a. In this case, coincides with that of the 
(a, 1,1)-[DK] model, but 

2 Xoo(l - Xoo) (a^ - (a^ + 2a - 1) Xoo) 



cr 



(a - (1 + a)xoo)^ 

This generalized form of the [MT] model in which a spreader becomes a stifler only after 
being involved in a random number of unsuccessful communications is further studied in 
Lebensztayn et al. jsj. 



Example 2.10. Other rumour models in the literature for which our results apply were 
proposed by Pearce l3| and Kawachi f^. The first one has the dynamics of the [DK] 



model, with the following interaction rules. A meeting between an ignorant and a spreader 
results in the ignorant turning into a spreader with probability p. When two spreaders 
interact, either both become stiflers with probability q2 or only one of them does so with 
probability q\ < 1 — ^2- Finally, the interaction of a spreader with a stifler results in two 
stiflers with probability r. Thus, this model is obtained by considering X = p, 6 = 1, 
01 = (I2/P, O2 = qi/ (2p) and 7 = r/p. 

One of the deterministic models studied in Kawachi |6] evolves similarly to Hayes' model, 
as explained in Example 12.71 The contacts are as in the [MT] model, but, when two spread- 
ers meet, both become stiflers with probability (3; the transformation of only one of them 
into a stifler is not possible. When a spreader encounters an ignorant, the first one trans- 
mits the rumour with probability a and the latter joins the spreaders with probability 6. 
Finally, in a meeting between a spreader and a stifler, the first individual becomes a stifler 
with probability 7. This model is obtained by the choice X = a, 6 = 9, 9i = 2 f3/a, 62 = 
and 7 = 7/0. 

3. Proofs 

The main method for proving Theorems 12.41 and 12.51 is, by means of a random time 
change, to deflne a new process {V^^\t)}t>o with the same transitions as {V^'^\t)}t>o, so 
that they terminate at the same point. This transformation is done in such a way that 
{V^^\t)}t>Q is a density dependent Markov chain to which we can apply Theorem 11.4.1 
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of Ethier and Kurtz |4] . The arguments are similar to those used in Kurtz et al. |2J for the 
coverage of a random walks system on the complete graph. 

3.1. Random time change. We define 

Jo 



poo 

T(^)(s) = inf{t : e(^)(t) >s},0<s< / F^^^ 

Jo 



u) du, 



and let V^^\t) = V^^\T^^\t)). The time-changed process {V"*^^^(t)}(>o has the same 
transitions as {V^^\t)}t>Q, hence if we define 

f(^) = inf{t:F(^)(t)=0}, 

we have that V^(^)(r*^^)) = \/(^)(f(^)). Furthermore, {V^'^\t)}t>o is a continuous-time 
Markov chain with initial state (A^, 0, 1) and transition rates given by 

transition rate 

4 = (-1,0,1) X6X, 

(3.1) £i = (-1,1,0) X{1-6)X, 

^2 = (0,0, -2) ^^1^-^' 
£3 = (0,0,-1) \e2iY -l) + \-f{N+l-X -Y). 
3.2. Deterministic limit of the time-changed process. We define for t > 0, 

^(^Ht) = ^^ = (x(^)(t),n(^)(t),^(^)(t)), 

and consider 

f3eo (x, u,y) = XSx, f3e^ {x,u,y) = X{1 - 6) x, 

y 

f3(,^{x,u,y) = \ 6i-, P(,^{x,u,y) = A 6*2 y + A 7 (1 - x - y). 

Note that the rates in (13.11) can be written as 

N 

so {v^^\t)}t>Q is a density dependent Markov chain with possible transitions in the set 

Now we use Theorem 11.2.1 of Ethier and Kurtz to conclude that the time-changed 
system converges almost surely as — )■ 00 (on a suitable probability space). The drift 
function is given by 

F{x,u,y) = ^ii(3t,Xx,u,y) = (-Ax, A(l - 5)x, A(7 + 5)x - \9y - A7), 



^ . X U Y\ ^ / 1 
iV'iV'iV +^liV 
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hence the hmiting deterministic system is governed by the following system of ordinary 
differential equations 

'x'(t) = -Ax(t), 
u'{t) = X{l-6)x{t), 
' y'{t) = \{j + 6)x{t)-\9yit)-\j, 
^x(O) = 1,m(0) = and y{0) = 0. 

The solution of this system is given by v(t) = {x(t),u(t),y{t)), where 

= e-^*, u{t) = (1 - - x{t)), y{t) = f{x{t)), 

with / replaced by fe if 6 equals or 1. According to Theorem 11.2.1 of Ethier and 
Kurtz Q, 

Lemma 3.1. We have that v^^\t) converges almost surely to v{t), uniformly on bounded 
time intervals. 

We prove that for each of the first two components of v^^\ the convergence is uniform 
on the whole line. 

Lemma 3.2. We have that x^^\t) converges almost surely to x{t), uniformly on M. The 
analogous assertion holds for u^^^ and u. 

Proof. We prove the first statement. Given any e > 0, we take to = ^o(^) such that 
x{t) < e/2 for all t > to- By the uniform convergence of x^^\t) to x{t) on the interval 
[0,to], there exists A^o = ^o(^) such that 

\x^^\t) - x{t)\ < e/2 for all > A^o and t G [0,to]. 

Therefore, for all N > Nq and t > to, 

^(iv)(i) < x^^\to) < x{to) + e/2 < £, 

whence \x^^\t) — x(t)| < e. □ 

3.3. Proofs of Theorems 12.41 and 12.51 Both theorems follow from Theorem 11.4.1 of 
Ethier and Kurtz [4.]. We adopt the notations used there, except for the Gaussian process 
V defined on p. 458, that we would rather denote by V = (V^, Vu, Vy). Here ip{x, u, y) = y, 
and 

Too = inf{t : y{t) < 0} = -- logXoo- 

Moreover, from (12.11) . 

(3.2) V^{v{t^)) ■ F{v{t^)) = y'{T^) = A(7 + 5)x^-\-i< 0. 

Let us explain the argument leading to the proof of the Law of Large Numbers. Although 
1/(0) = 0, we have that y'if)) > 0. This and (13. 2p imply that yijao — s) > and y{Too+£) < 
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for < e < Too- The almost sure convergence of y^^^ to y uniformly on bounded intervals 
yields that 

(3.3) lim f = Too almost surely. 

Thus, recalling that V^^\t^^'^) = \/(^)(f(^)), we obtain Theorem 12.41 from formula (13. 3 p 
and Lemma 13.21 

With respect to the Central Limit Theorem, we get from Theorem 11.4.1 of Ethier and 
Kurtz [J] that y/N (x^^-'(f^^^) — Xoo, S*^^-*(f (^^) — Mqo) converges in distribution as — ?■ oo 
to 

(3.4) (V,(r^) - A Vy{T^), Vu{r^) + A{l-6) V,(r«,)) , 



where A is the constant defined in Theorem 12. 51 The asymptotic distribution is a mean zero 
bivariate normal distribution, so it remains to explain how formula (12. 4p for the covariance 
matrix S is obtained. 

For this, we summarize the steps taken to compute A = Cov(V(roo), V(roo)), a task 
that can be better carried out with mathematical software. First, we calculate the matrix 
of partial derivatives of the drift function F and the matrix G, as defined in Ethier and 
Kurtz 0, p. 458]. They are given by 

/ -A \ 

dF{x,u,y) = I X{l-6) and 

\A(7 + 5) -xej 

(\x — A(l — 5)x —X5x 

-X{l-6)x X{l-6)x 
-X6x X{6 - -f)x + X{k - e + 2-^)y + X-f ^ 

Next, we obtain the solution $ of the matrix equation 

^ $(t, s) = dF{x{t),u{t),y{t)) $(t, s), $(s, s) = h, 
where Is is the 3x3 identity matrix. Then, we compute 

(3.5) Cov(V(t), V{t)) = f <l>(t, s) G{x{s),u{s), y{s)) [$(t, s)f ds. 



We emphasize that the computation of A must be separated into four cases: 9 G (0, 1) \ 
{1/2}, 9 = 1/2, 9 = and 9=1. The value 9 = 1/2 must be considered separately from 
the interval (0, 1) owing to the appearance of the integral e'*'^^^"^^'' ds in the element (3, 3) 
of the matrix given in formula (13. 5p . This also explains why the constant D in Theorem l2.5l 
is defined differently for 9 = 1/2. 

The final step to obtain A is to set t = Too in the formula obtained from (13. 5p . This is 
accomplished by making suitable substitutions in this formula according to the value of 9; 
for instance, for 9 G (0, 1) \ {1/2} we replace e~^^ and e~^^^ respectively by 

7 + 06' 
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The resulting formula (valid for any 6) is 



(3.6) 



A 



v 



-(1 - 5) (1 - Xoo) Xc 







-{1 - 5) (1 - x^) (1 - 5) (1- Xoo) ((1-5) + 5) -B\. 
-B D 



Using (13.41) . (13. 6 p and the well-known properties of the variance and covariance, we get 
formula (I23D- 
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